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ABSTRACT 

Aims. Here we report evidence of a large solar filament eruption on 2013, September 29. This smooth eruption, which passed without 
any previous fiare, formed after a two-ribbon fiare and a coronal mass ejection towards Earth. The coronal mass ejection generated 
a moderate geomagnetic storm on 2013, October 2 with very serious localized effects. The whole event passed unnoticed to flare- 
warning systems. 

Methods. We have conducted multi-wavelength analyses of the Solar Dynamics Observatory through Atmospheric Imaging Assembly 
(AIA) and Helioseismic and Magnetic Imager (HMI) data. The AIA data on 304, 193, 211, and 94 A sample the transition region and 
the corona, respectively, while HMI provides photospheric magnetograms, continuum, and linear polarization data, in addition to the 
fully inverted data provided by HMI. 

Results. This flux emergence happened very close to a filament barb that was very active in mass motion, as seen in 304 A images. 
The observed flux emergence exhibited hectogauss values. The flux emergence extent appeared just beneath the filament, and the 
filament rose during the following hours. The emergence acquired a size of 33" in ~12 h, about ~ 0.16 km s“k The rate of signed 
magnetic flux is around 2x10^^ Mx min"^ for each polarity. We have also studied the eruption speed, size, and dynamics. The mean 
velocity of the rising filament during the ~ 40 min previous to the flare is 115+5 km s~\ and the subsequent acceleration in this period 
is 0.049+0.001 km s“^. 

Conclusions. We have observed a supergranular-sized emergence close to a large filament in the boundary of the active region 
NOAA11850. Filament dynamics and magnetogram results suggest that the magnetic flux emergence takes place in the photospheric 
level below the filament. Reconnection occurs underneath the filament between the dipped lines that support the filament and the 
supergranular emergence. The very smooth ascent is probably caused by this emergence and torus instability may play a fundamental 
role, which is helped by the emergence. 
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1. Introduction 

On 2013, September 29 a large solar filament rose smoothly, 
without any previous flare, vigorous dynamics, or significant 
brightening arcade. This filament lifted and afterwards created a 
two-ribbon flare, and its X-ray flux reached the category Cl.2 by 
GOES 1-8 A^ The eruption led to an Earth-directed halo coro¬ 
nal mass ejection (CME). The interplanetary plasma interacted 
with the Earth’s magnetic field creating a moderate geomagnetic 
storm with effects on the ground (Dst geomagnetic index reached 
-75 nT, and storm classification follows Gonzalez et al. 1994) 
on 2013 October 2^. However, since the interplanetary medium 
and terrestrial effects are beyond the scope of this paper, they will 
be treated elsewhere. This eruption went unnoticed for flare early 
alerts and predictions and it produced a geomagnetic storm. This 
case is a very pertinent example, which can be used to obtain in¬ 
put for simulations and for solar and space weather predictions 
and warnings. 

In this paper, we focus on the solar source and trigger of this 
filament eruption, which seems to be a supergranular-sized mag- 


^ http: //WWW . Imsal.com/solarsoft/last_events_201310(S)l_ 
2320/index.html 

^ Dst data available in http: //wdc. kugi. kyoto-u. ac. jp/dst_ 
provisional/201310/index.html 


netic flux emergence occurrence underneath the filament. We in¬ 
vestigate the reasons for this kind of an undisturbed rise using 
Solar Dynamics Observatory (SDO) data by means of the Atmo¬ 
spheric Imaging Assembly (AIA) at various EUV wavelengths 
(AIA/SDO 304 A, 193 A, 211 A and 94 A). Also, to explain the 
event we obtained a close-up using Helioseismic and Magnetic 
Imager (HMI) magnetograms, continuum, and linear polariza¬ 
tion images. We also took advantage of the fully inverted HMI 
data in this work. 

In the review of Linker et al. (2003) it is stated that only 20% 
of the CMEs are associated with a large flare. Elareless CMEs 
have been reported by Song et al. (2013). Gosain et al. (2012) 
state that most of the filament eruptions related to CMEs pass 
undetected because of their low emission and, consequently, are 
not associated with flares. Therefore, trusting only in flare occur¬ 
rence may be misleading and the region can actually be CME- 
productive. 

Some very recent and comprehensive summaries of the vari¬ 
ous models that can trigger the eruption of filaments and CMEs 
are Aulanier et al. (2010); Schmieder et al. (2013); Aulanier 
(2014); Parent! (2014) and references therein; the "magnetic 
breakout model" (Antiochos 1998; Antiochos et al. 1999); a 
highly twisted flux rope, as in the "torus instability" (Torok & 
Kliem 2003, 2005; Kliem & Torok 2006); different instabilities. 
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e.g. MHD (Klimchuk & Sturrock 1989; Sturrock 1989); catas¬ 
trophic loss of equilibrium (Forbes & Isenberg 1991); approach¬ 
ing polarities (Forbes & Priest 1995); moving features and mag¬ 
netic cancellation (van Ballegooijen & Martens 1989; Forbes & 
Isenberg 1991); and the "tether cutting" model, implying recon¬ 
nection in the arcades below the flux rope (Moore & Roume- 
liotis 1992; Moore et al. 2001). Considering the magnetic flux 
emergence, filament eruptions due to flux emergence in the sur¬ 
roundings have been analysed in simulations by Chen & Shibata 
(2000). The magnetic configuration of active and quiescent fil¬ 
aments and their CME productivity have been studied in Feyn¬ 
man & Martin (1995). Slow and fast eruptive phases have been 
studied by Sterling & Moore (2005); Sterling et al. (2007). We 
consider the observational features of these models to choose the 
most plausible model for this phenomenon. 

More particularly, different works are revised in the review 
by Chen (201 1): Zhang et al. (2008) studied events of flux emer¬ 
gence and cancellation 12 h prior to a CME initiation, finding 
that 91% are related. Wang & Sheeley (1999) also found filament 
eruptions close to flux emergence regions. Jing et al. (2004) de¬ 
tected that more than half of the eruptive filaments were related 
to flux emergence. Zhang et al. (2001) related relatively small 
flux cancellation areas with the filament eruption of the Bastille 
Day event. 

However, many of these works were qualitative, only with 
the presence and closeness of the flux emergence to relate with 
other magnitudes. Here we aim to describe the filament quanti¬ 
tatively and to assess its instabilities, and we describe how a flux 
emergence can influence a filament lift-off. 

This paper is organized as follows: first, the observations are 
described in Section 2, describing the flux emergence and fila¬ 
ment dynamics, along with the surrounding coronal holes. All 
calculations are included in Section 2. The discussion and con¬ 
clusions are presented in Sections 3 and 4, respectively. 

2. Observations of the emerging region and 
fiiament iift-off 

We used data provided by the space solar facility Solar Dy¬ 
namics Observatory (SDO, Pesnell et al. 2012). Data were ac¬ 
quired with the Atmospheric Imaging Assembly (AIA, Lemen 
et al. 2012) and the Helioseismic and Magnetic Imager (HMI, 
Scherrer et al. 2012). The AIA multi-wavelength data comprise 
observations in 304 and 193 A, which sample the high chro¬ 
mosphere, transition region, and corona, respectively. We also 
included AIA 211 A data, which are useful as an intermedi¬ 
ate wavelength measurement since the filament opacity, com¬ 
pared to 304 A, is somewhat reduced. The AIA 94 A data also 
probe the inner corona. The HMI data sample the Fe i line at 
6173 A to generate line-of-sight (hereafter, LOS) photospheric 
magnetograms. This line samples a photospheric height around 
~100 km (Fleck et al. 2011). The data cadence for the general 
study is chosen to be 12 min from September 29 to October 1 on 
both AIA and HMI instruments. In other particular cases of the 
study, the cadence is reduced to 1 min or increased up to 180 min 
(3 h). Images are co-spatial and quasi-synchronized (within sec¬ 
onds). Pixel sampling for both instruments is 0"6. 

All images are reduced through the SolarSoftWare (SSW) 
read_sdo.pro and aia_j?rep.pro routines. Next, frames are cor¬ 
rected from rotation and aligned. 

Figure 1 shows context images. The left panel indicates the 
whole solar hemisphere, showing the filament and two surround¬ 
ing coronal holes in AIA 193 A. The right panel shows the fila¬ 


ment in AIA 304 A as a dark, inverse-S-like, elongated structure 
located between two opposite polarities of NOAA ARl 1850 and 
the following active region. The filament was visible close to 
the trailing boundary of AR 11850, with its southern part close 
to its positive facular polarity, and the spine located at the po¬ 
larity inversion line (PIL), with a facular negative polarity on 
the left and another positive on the right. As seen in AIA 304 
and 193 A, in close to five days prior to the eruption, the fila¬ 
ment developed from a short structure (Sept 25, 21:00 UT) to a 
inverse-S-shaped elongated in just one day. This fully developed 
filament evolved in a dark area in 193 A, which is an open held 
region. The length was about -0.407?© on Sept 25 to -0.637?© 
two days later. This filament showed very active mass motion 
along the spine and barb. The southern part is more elevated 
than the northern part at this stage. The heights look similar at 
Sept 28, 10:00 UT. This situation seems to reverse on Sept 28 
19:00 UT. The magnetic flux emergence started around Sept 29 
at 00:00 UT and the appearance of pushing the filament from 
below is conspicuous. This emergence displayed small pores in 
continuum. There was a positive polarity where one of the barbs 
was located, and this region seems to be attaching the filament 
to the surface. The apex of the filament rose, placed somewhat 
northern from the barb, around 21:00 UT. Eventually the fila¬ 
ment developed a two-ribbon flare on Sept 29 at 21:43 UT, which 
led to a maximum flux of X-ray of Cl.2 in GOES. The southern 
leg was detached early in the flare, and the northern leg produced 
a shower of coronal rain on Sept 29, at 23:29 UT, lasting around 
one hour. Thread-like features of the filament feet appeared to 
fall down. 

The subsequent CME appeared in LASCO at 22:24 UT, 
with a plane-of-the-sky speed (POS) of - 600 km s“^ . 
The halo CME is reported in the CACTUS catalog 
(http://side.oma.be/cactus/catalog/LASC0/2_5_0/ 
qkl/2013/09/latestCMEs.html, Robbrecht & Berghmans 
2004; Robbrecht et al. 2009). However, the reported speed in the 
LASCO catalog (http : //edaw. gsfc. nasa. gov/CME_li st/) 
is 1179 km s“^ This discrepancy is due to the external CME 
shell, which is faster than the filament. 

2.1. Magnetic flux emergence 

We selected AIA 304 A images and HMI magnetograms with a 
cadence of 12 min, from Sept 29 00:00 UT to Oct 1 10:00 UT. 
We analysed AIA 304 A close-up images with overlaid HMI 
magnetograms of the aforementioned filament barb (located on 
- 18°N, 15°W in the left panel of Fig. 1), which appears to be 
the location where the photospheric magnetic flux emergence 
happens. The panels in Fig. 2 show the evolution, displaying se¬ 
lected frames^. This barb seems to be close to a positive po¬ 
larity, on Sept 29 at 00:00 UT. A large magnetic supergranular 
extent appeared, as seen in the HMI magnetograms. This large 
supergranular emergence exhibited mixed polarities at first, and 
later on the positive and negative polarities drifted, rearrang¬ 
ing as ordinary network elements. Three hours later, the main 
negative polarity started developing a pore. The positive polari¬ 
ties also started developing small pores in the following hours. 
Six hours later, this emergence extent, consisting of positive po¬ 
larities (green) on the left and negative polarities (blue) on the 
right, moved apart up to - 25" (panel 1 in Fig. 2), and reached a 
maximum size of - 33" around 13:00 UT. At this moment, the 
negative polarity peaked -1000 G, while the positive was about 

^ The temporal evolution shown in Fig. 2 is available in the on-line 
edition. 
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Fig. 1. Context figures of the filament. The grid spaces mark 15°. Left: Image in AIA 193 A, on 2013 Sept 29, at 00:00 UT. Image is clipped 
from 50 to 5000 DN. The filament is located on 20°N 15°W, with two coronal holes surrounding it on coordinates 40°N 20°W (CHI) and 0°N 
7°W (CH2). Right: Context image showing a derotated AIA 304 A image displaying the filament, and the positive polarity (in green) and negative 
(blue). Different colour shades are explained in caption of Fig. 2. The bipolar flux emergence occurrence is located at [250", 200"], enclosed in 
a white square. A movie available in the online data base shows the temporal evolution of the whole visible hemisphere from 21:00 to 22:09 UT 
in the AIA 304 A channel. In the movie, the FOV of the right panel of this figure is marked by a black rectangle, and the emergence by a white 
square. The green and blue contours have the same meaning as in the right panel. 


800 G. Therefore, the magnetic field of both polarities increased 
from around +(-)200 to +(-)900 G in only six hours. A negative 
polarity on Sept 30, 11:24 UT, emerged and drifted to the right 
side with other negative polarity elements. These magnetic re¬ 
arrangements are visible for the whole period and the loops in 
AIA 304 A mark the polarity correspondence; however, no pro¬ 
nounced change in polarity is noticeable in HMI magnetograms, 
neither shear in AIA 304 nor 193 A data. At the end of the se¬ 
quence, the polarities diminished their flux and approached each 
other. The magnetic emergence seems to push the barb of the 
filament from below during the whole event. 

In the second row of panels in Fig. 2, we plot AIA 211 A im¬ 
ages in mainly two lines: Fe xiv, corresponding to logT ~ 6.30; 
and Fe xiii, corresponding to logT ~ 6.25 (O’Dwyer et al. 2010). 
This wavelength happened to become very adequate to observe 
flux emergence patches under filaments, since the S/N is good, 
samples high temperatures, and the opacity of the filament be¬ 
comes reduced but not totally suppressed. 

We also present evidence of an EUV counterpart of the flux 
emergence in 94 A in the third row of panels in Fig. 2. These 
images sample a temperature range in the corona (O’ Dwyer et al. 
2010) with two lines: the inner quiet corona, sampled by Fe x, 
corresponding to logT ~ 6.05, and the flaring range, sampled 
by Fe xviii, corresponding to logT ~ 6.85. The region started to 
increase its intensity from 6 DN to 30 DN. In O’Dwyer et al. 
(2010), the quiet sun intensity is synthetically defined as 6 DN. 
Therefore, the intensity of this emerging region is five times that 
of the quiet sun, and 10% of the maximum intensity in the flare 
of the full disk in this period. 

This observational fact may suggest rapid connectivity be¬ 
tween the photospheric longitudinal and transverse fields and 
the corona, since the signal starts to increase from 01:00 UT on¬ 
wards, as shown in Fig. 2 and in the movie added as on-line 
material. 


After the flare, we can see that opposite magnetic polari¬ 
ties rearrange and sometimes approach and disappear, display¬ 
ing a brightening between them as a cancellation signature. In 
the post-flare panels, the last one is a very clear case. 

Since we want to measure the magnetic flux emergence and 
check whether magnetic cancellation takes place, we compute 
the magnetic flux and magnetic flux density to this small area 
in HMI FOS-magnetograms. We set a threshold of +(-)90 G, 
that is, lOcr for 45-s magnetograms (Fiu et al. 2012). This 
threshold is high enough to account for the flux computations in 
the inter-network and network, avoiding noise effects. Under all 
considerations, the Ailing factor is equal to 1. 

We calculated the signed magnetic (positive and negative) 
flux separately, as in the left panel of Fig. 3. The magnetic flux 
has been corrected from heliocentric angle (yu). The vertical bars 
mark the flare onset. The starting time is Sept 29 at 00:00 UT. 
The positive flux reaches a maximum value at Ar=600 min, but 
there is a plateau from A^=600 min to A^=960 min (from 10- 
16h UT). Both positive and negative fluxes start to decrease 
around Sept 29 at 12:00 UT. The maximum unsigned flux is 
~3.4xl0^^ Mx, and the unsigned flux rate, 4.29+0.01 xlO^^ Mx 
min“^. Calculating the emergence flux rate (up to the middle of 
the plateau) yields + 2.23+0.07 xlO^^ Mx min“^ Considering 
the first maximum, the rate is +2.60+0.08 xlO^^ Mx min“^ We 
performed the same operation for the negative polarity, obtain¬ 
ing -2.52+0.07 xlO^^ Mx min“^ For the negative polarity, we 
can see a clear decay in both magnetic flux and flux density, as 
shown in Fig. 3. Calculating the decay rate of this negative area 
yields 5.2+1.2 xlO^^ Mx min“\ which is much slower than the 
emergence rate. 

The signed magnetic flux densities, shown in the right panel 
of Fig. 3, are both around ~+(-)250 G and follow the trend as 
the signed magnetic flux. There are two peaks on the positive 
magnetic flux density some time after the flare, probably caused 
by recurrent magnetic intensification in the area (see Section 3). 
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To rule out any possible contribution from other areas to the 
flux emergence in the region, we have also analysed three dif¬ 
ferent areas where magnetic emergence is detected (close to the 
main sunspot of the active region [700", 90"]. This area was ac¬ 
tually flare-productive on Sept 29 at 05:11 UT, with the onset of 
a Cl.6 flare. In the active region, the positive polarity decreased 
from 2x10^^ Mx about 25%, and the negative polarity increased 
from -1.5x10^^ Mx by 16%, becoming both fluxes equal 2.5 h 
before the flare. Furthermore, we studied the areas around coor¬ 
dinates [400", 150"] and [400", 30"], corresponding to super- 
granular voids, where magnetic flux emergence is also observed. 
For these regions and the threshold set, the flux emergence rate is 
negligible. We computed the total magnetic flux of the whole AR 
as well, from September 24 00:00 UT to Sept 29 00:00 UT. The 
positive flux started at 0.9x10^^ Mx, peaking on Sept 26 with 
1.3x10^^ Mx. The negative flux started at -1.1x10^^ Mx, peak¬ 
ing at the same time. At the time of the flare, the fluxes were 
similar to Sept 24, keeping the decreasing trend. The rate of the 
total flux emergence respect to the AR is around 1%. 

We also measured the expansion velocity of the emergence, 
with two main opposite polarities, from Sept 29 01:00 UT to 
13:00 UT with different methods. The first method is the mag¬ 
netic centroids, as detailed in Balmaceda et al. (2010), but the 
distance was computed between the positive and negative po¬ 
larities above the threshold of +(-)90 G. With this method, the 
expansion velocity is 0.165+0.006 km s“^ The second method 
we used is visually measuring the distance between the two 
main polarities (as a diameter and obtaining the radius), getting 
0.100+0.003 km s“^ The third procedure is the one described 
in Palacios et al. (2012), measuring the area and assuming it as 
circular, subsequently calculating the radius from it. This yields 
a velocity of 0.157+0.004 km s“^ 


2.1.1. Linear polarization and transverse field of the flux 
emergence 

SDO/HMI retrieves a full-disk imaging of Stokes parameters, in 
the form of HMI full-disk Stokes spectrograms (HMI.S _720s 
data series). We used Stokes U and Q images to create linear 
polarization images. These images are created with six points 
(line sampling of about 69 mA from the centre line; Schou et al. 
2012; Centeno et al. 2011) in the Fe i line at 6173 A. Using the 
standard SSW routine for HMI/SP reduction, images are accu¬ 
mulated and processed. Profiles are extracted from the images 
and normalized by the red-most wavelength point in Stokes I 
(J^). For each frame, we have computed the linear polarization 
using the following formula (e.g. Solanki et al. 2010; Martmez 
Pillet etal. 2011): 



z 



( 1 ) 


In this stage, we also used also HMI inverted data by the 
Milne-Eddington inversion code VFISV (Borrero et al. 2011; 
Centeno et al. 2014). These data consist of the magnetic held, 
azimuth, and inclination to compute and plot the transverse held 
through B • siniy) with y the inclination with respect to the LOS, 
as described in Hoeksema et al. (2014). For a general comparison 
between the LOS-magnetograms and the LOS-magnetic held 
obtained by weak held approximation through the Stokes pro¬ 
files in active regions, see Palacios et al. (2014). 

We studied the transverse held magnetic flux density. This 
includes a threshold of 180 G, which is double the longitudinal 


flux. The magnetic flux density is included in Fig. 4. In the left 
and central panels, we show both longitudinal (green and blue) 
and transverse held (shades of orange), with a 94 A background, 
in two different instants: 03:48 and 12:00 UT. The rightmost 
panel shows the linear polarization at 03:48 UT plotted in loga¬ 
rithmic colour scale. The maximum linear polarization is 1.5%. 
The structure of the transverse held in the blob is not uniform, 
but creates bridges between the main positive and negative polar¬ 
ities. In the first four hours, this situation is conspicuous. After 
that, the transverse held patches remain in the main polarities 
because some inclination of the mainly vertical magnetic held. 
Some showers of cosmic rays appear at some moments, altering 
the most sensitive part of the polarization images, i.e. the linear 
polarization, translated into transverse magnetic held and mak¬ 
ing it look noisy. 


2.2. Filament rising 

We calculated the rising speed on the apex of the filament on 
de-rotated AIA 304 A images, also provided in the on-line ma¬ 
terial, and whose detailed reference frame is shown in the right 
panel of Fig. 1. The cadence of this series is 1 min from 21:00 
to 22:00 UT. We examined the previous frames from 20:00- 
21:00 UT, but the rising is negligible at that period. As shown 
in Fig. 5, the filament kept a constant height during 22 min, and 
the lift-off took around 38 min.The lift-off of the filament started 
much before the flare onset at 21:43 UT. The height was esti¬ 
mated on de-rotated images, considering a fiducial point [340", 
340"] that lies in the line of the axis defined by the two-ribbon 
flare middle point and the footpoints. The distance between foot- 
points is around 0.8 Ro.The apex height was measured, albeit 
the southern part of the filament was more active on rising the 
first 20 min. These measurements are indicated with triangles in 
Fig. 5. The visual error for the uncorrected height is considered 
as 3 pixels of a rebinned image of 1024 x 1024 pixels, translated 
as 5220 km (not shown). We corrected from projection of the 
heliocentric angle, i.e. dividing by the heliocentric angle sine, 
considering the eruption follows the local surface normal direc¬ 
tion, as shown in the plots of Fig. 5. The mean velocity from the 
flare to the end at 22:00 UT is 188+5 km s~\ and the subsequent 
acceleration in this period is 0.065+0.007 km s“^. The maximum 
velocity is 312 km s“^ (Considering the 38-min period that the 
rising lasts, the velocity is 115+5 km s“^ and acceleration is 
0.049+0.001 kms-2). 

For this event, we fit the ascent height to different func¬ 
tions, since this can shed light on the physical mechanism 
that triggered the eruption. We tried exponential, log - log, and 
parabolic function. Then we evaluated them with the merit func¬ 
tion (unreduced). The exponential fit for the time-height 
(uncorrected) yields ;^^=1.7, while the log - log fitting yields 
;^^=5.9. In the case of projection-corrected height, the fitting is 
even better: ;^^=0.9, while the log - log fitting yields ;^^=3.4. 
Parabolic fittings get very large 

We also completed the study with SOHO (Domingo et al. 
1995) LASCO C2 and C3 level 1 data (Brueckner et al. 1995). 
After centring and rotating, we performed image base differ¬ 
ences and then, applied a Sobel operator for enhancement. We 
measured the apex in the filament images. The plate scale for C2 
is 11.9"pix“^ and for C3 is 56"pix“^ We show the whole time- 
height scale for AIA304, C2 and C3 in Fig 5 (bottom). The;^^ 
for the whole sequence is 3.80. 
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Fig. 2. Top row. frame sequence of AIA 304 A images of a filament’s barb (darker areas) with LOS-magnetic field superimposed. Different shades 
of green (from lighter to darker) contours mark different levels, as 90, 200, 500, 600, and 700 G, while a gamut of blue contours mark -90, -200, 
-500, -700, -900 G for the negative polarities. The frame size is 100"x 100". Sequence timing aims to show the magnetic flux emergence (panel 
1), pre-fiare (panel 2), flare (panel 3), and decay (panel 4), from left to right. Bottom row: same temporal sequence in AIA 211 and 94 A. Movie 
issued as on-line material. 




Fig. 3. Left: computed magnetic flux on positive patches (solid line), on negative patches (dashed line), signed magnetic flux (dotted-dashed line), 
unsigned magnetic flux (dotted line). A vertical line marks the time of the flare at 21:43 UT. Right: mean values of the magnetic flux density of 
the positive (solid) and negative (dashed) regions. 
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Fig. 4. Left and centre: AIA 94 A image at different instants with LOS-magnetic field superimposed, with the same contours as Fig. 2. Orange and 
red contours mark levels of 180, 200, 250, 400, 500 G in transverse field. Right: linear polarization map at 03:48 UT, plot in logarithmic colour 
scale. 
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Fig. 5. Left: height versus time of the filament apex. Two plots are shown: triangles represent the uncorrected height data; and the solid line 
indicates the height after correcting of heliocentric angle projection. Error bars are computed considering the visual error on height and angle 
correction. Time starts from 21:00 UT and covers 1 h. In the last 38 min, the height starts increasing. A vertical bar marks the flare onset, at 
21:43 UT. Right: height-time diagram, including LASCO C2 (diamonds) and C3 (asterisks) data. 


2.3. Potential field modelling and critical decay index for the 
filament region and emergence 


carello et al. (2014, and references therein), we estimated this 
magnitude as 


We used the package Potential Field Source Surface available 
in SSW PFSS (Schrijver & De Rosa 2003) to assess the critical 
height of the filament through the decay index of the external 
field n, since we aim to know whether the filament was already 
unstable before the emergence or the emergence made the fila¬ 
ment unstable. The underlying model for the decay index n is the 
torus instability. The input of the PFSS model is a synoptic mag¬ 
netogram. We used the highest resolution input data, the avail¬ 
able daily quasi-synoptic HMI magnetograms, obtained at 12:00 
UT, on Sept 25, 26, 28, and 29 (27 not available). These synop¬ 
tic maps were not derotated, so care had to be taken to identify 
the areas. We set the PFSS computation grid as 800 points in so¬ 
lar longitude, 400 in latitude, and 66 points in height, covering 
from IRq to 2Rq. Therefore, the spatial sampling of the grid is 
0.45° per grid point in longitude and latitude, and around 10600 
km per grid point on height. We included a context image of the 
PFSS loop set for this region on Sept 29, 18:00 UT (top panel of 
Fig. 6). The height rendered is in agreement with the estimated 
filament height. 

Then, with the decay index of the external field n, used re¬ 
cently in Kliem & Torok (2006); Torok & Kliem (2007); Zuc- 


n = -R 


dlnjBe,) 

OR 


( 2 ) 


We generated datacubes of n for Sept 25, 26, 28, and 29. 
We selected a box of 10x65x66 grid points, equivalent to 4° 5 
x29°25 xIRq (solar longitude, latitude, radial height). The Z 
axis actually goes from ~ 20 Mm to 700 Mm, since we removed 
the points on the very bottom. This final volume of ~55 Mm x 
357 Mm X 680 Mm includes the filament and flux emergence. 
The Y axis covers approximately the filament’s length, but con¬ 
sidering its shape, it cannot be considered as wholly lying in the 
solar longitude selected plane. We clipped the n values from 1.2 
to 2.0. We cut this box along and across the exact location where 
the flux emergence happens, as shown in Fig. 6. 

The decay index n is represented in shades of purple-blue 
on the volume where it is torus-unstable (n becomes a critical 
factor for instability, ricrit, when values ranges from 1.3 to 1.5, as 
in Zuccarello et al. 2014). Considering the low impulsivity of the 
ejection, this threshold is also in agreement with Torok & Kliem 
(2007). 

During the previous days, some network points are rooted 
in the area, but the most important variation are from Sept 28 
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12:00 UT to Sept 29 12:00 UT. The height of the filament may 
vary and is only marked for Sept 29, with a hollow circle the 
projection-uncorrected height, and the corrected height with a 
cross. On Sept 28, n increased on the emergence area. Black 
areas are those where ^ <1.2, and these areas decreased from 25 
onwards. On Sept 26, a large network area also made n increase, 
but it was not located in same coordinates of the emergence. 

3. Discussion 

The emergence of this relatively small-scale region follows 
Hale-Nicholson’s law of polarities (the same leading and trail¬ 
ing polarity in the same hemisphere. Hale et al. 1919); in this 
case, negative is leading, and positive is trailing, as the leading 
ARl 1850. However, in the right panel of Fig. 1, we notice a dif¬ 
ferent magnetic configuration. The filament is located between 
the negative polarity of the next active region (to the left) and 
the positive polarity from ARl 1850 between where we can lo¬ 
cate the polarity inversion line (PIT). Therefore, the emergence 
exhibits a normal configuration, but inverse to the filament con¬ 
figuration, which can help the magnetic instability. The situation 
is similar to the simulations that originate eruptions in Chen & 
Shibata (2000), as an opposite polarity configuration emerges 
beneath or at the side of a filament. Observing in different coro¬ 
nal wavelengths, as 193 A and 94 A, there is no evidence of a 
large magnetic canopy over the filament in the form of bright 
loops. LASCO images also reveal a leading edge and the fila¬ 
ment inside this edge. These images seem to be in agreement 
with the standard model of a twisted flux rope surrounded by a 
less twisted longitudinal magnetic field, but without including 
a large-scale magnetic arcade. The inverse-S shape and barbs 
would define it as dextral chirality, and the longitudinal field of 
the filament would then point to the solar north direction. The 
twisted field lines would exhibit negative helicity. 

We have to note that this filament configuration (dextral) fol¬ 
lows the hemispheric pattern; however, exceptions to this rule are 
reported (e.g. see Mackay et al. 2014, and references therein). 
Furthermore, under a full polarimetry analysis, the chirality def¬ 
inition may change, from a topological orientation of features 
to a real estimation of the magnetic field direction and sense 
(Hanaoka & Sakurai 2014). 

The barb is very close to the positive polarity from before 
the supergranular flux emergence, and seems fundamental to the 
mass dynamics of the filament. Since the barb is very close to 
the emergence, it may be the point where the instability is trans¬ 
mitted. 

The rate of flux emergence, which is about 2x10^^ Mx 
min“\ is at least one order of magnitude larger than those re¬ 
ported by Palacios et al. (2012). However, the expansion veloc¬ 
ity of the flux emergence is much slower, at 0.16 km s“^ , than 
in Palacios et al., which is approximately 0.6 km s“^ .In this 
case, the higher magnetic fields may slow down the velocity of 
the upfiowing plasma. The expansion velocity calculated with 
centroids is similar. 

Expecting dramatic changes in photospheric magnetic flux 
as large reconnection features, we only found the bright blob in 
94 A, being conspicuously bright in 211 A. The magnetic flux 
emergence lasts around 12 hours and smoothly decays, proba¬ 
bly due to the reconnection that started a bit after the flux emer¬ 
gence detection. We included spectropolarimetric data to clar¬ 
ify the magnetic structure by adding the transverse component 
during the whole emergence. We may expect a larger transverse 
field structure, but that may be smeared out because of the av¬ 
eraging over 12 minutes, or some circular to linear cross-talk. 


Combined images between AIA 94 A images and longitudinal 
magnetograms suggest the magnetic link between the emergence 
in the photosphere and the coronal counterpart, in the form of a 
blob, brighter than the surroundings. This coronal counterpart is 
also visible in 304 and 193 A, but somewhat occulted by the barb 
due to its opacity. The blob is best seen in 211 A since the emer¬ 
gence patch is better seen through the filament barb, as the latter 
is less opaque. 

Compared to simulations of Kusano et al. (2012) and Tori- 
umi et al. (2013), this emergence region resembles an OP-type 
configuration with azimuth angle (pe ~ 160°, and shear ~ 
60° (0e is the angle between the PIT and the orientation of the 
bipolar region. With this shear and according to these simula¬ 
tions, the injected kinetic energy might be large enough to pro¬ 
voke an eruption. Bamba et al. (2013) states that OP-type may 
influence the destabilization of a flux rope via reconnection. 

We followed the evolution of the filament for about six 
days, and the phenomenon known as the "sliding door effect" 
(Okamoto et al. 2008, 2009; Kuckein et al. 2012) is not evi¬ 
dent, since the distance between the positive and negative facular 
boundaries remains the same, -220". The filament was already 
there at the beginning of the study, and the part that grew from 
there does not exhibit any distinctive pattern of the sliding door 
effect. 

Development of pores have been reported by Kuckein et al. 
(2012), and this is considered tracers of the filament emergence. 
However, Vargas Dominguez et al. (2012) had studied the event 
of Okamoto et al. (2008, 2009) and considered these photo- 
spheric signatures as necessary but not sufficient pieces of ev¬ 
idence for supporting the hypothesis of emergence. We also find 
very small appearance of pores. However, it is probably more 
related to the intensification of the magnetic field (Grossmann- 
Doerth et al. 1998; Bellot Rubio et al. 2001; Danilovic et al. 
2010) rather than characteristics of another flux rope emerging, 
since the magnetic field density increases from -200 to -1000 G 
. In the right panel of Fig. 3, other peaks appear in the positive 
magnetic field density, probably caused by recurrent magnetic 
field intensification in the region. 

Another point is that the CME speed was reported as 600 
km s“^. In the case of CME velocities, we always should have in 
mind the solar escape velocity, that is 617 kms“l Any difference 
from this velocity can account for the projection angle or the 
kinematic energy injection to the CME. 

Furthermore, we checked that CHI becomes darker visually, 
which might be an indicator of more open field lines there. Some 
of these features will be treated in another paper. 

Regarding the filament eruption, the current instability mod¬ 
els might explain or rule out the mechanism that causes the 
eruption. In the review of, e.g. Aulanier (2014), different mod¬ 
els are presented. The breakout model (Antiochos 1998; Antio- 
chos et al. 1999) depends on significant shearing, and involving 
reconnection above the flux rope (Schmieder et al. 2013). Re¬ 
connection is strongly enclosed in some of these models (Zuc- 
carello et al. 2013). Loss of equilibrium models (Forbes & Isen- 
berg 1991; Forbes & Priest 1995) may explain the slow fila¬ 
ment rising (Schmieder et al. 2013). These models rely on photo- 
spheric motions, polarity approaching, or vortical motion (Amari 
et al. 2003; Aulanier et al. 2010). Another model to consider is 
"tether-weakening", with observational examples presented by 
Sterling & Moore (2005); Sterling et al. (2007). This weakening 
is usually due to emergence off the filament channel (Moore & 
Roumeliotis 1992). 
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Fig. 6. Top left: PFSS magnetic rendering corresponding to the filament region rendered by PFSS for Sept 29 data. A low-resolution magnetogram 
was used, with the same FOV as Fig. 1 (right). Top right: transverse slice of the decay index n across the volume that encloses the filament 
and emergence. The estimated height of the filament for Sept 29 is marked with a black cross (corrected), or small circumference (uncorrected). 
Bottom: longitudinal slice of the decay index n along the volume that encloses the filament and emergence, and same symbols for the filament 
height. This figure belongs to subsection 2.3. 


In addition to the observational signatures of the models, the 
height-time fitting may provide hints about the physical mech¬ 
anism. In Schrijver et al. (2008), the expansion is fit with dif¬ 
ferent functional forms (polynomial, exponential, power-law...) 
and evaluated with the merit function Time-height exponen¬ 
tial fitting is also used in very rapid CMEs, like in the case of 
Gallagher et al. (2003). Gosain et al. (2012) also uses fitting to 
exponential. The time-height fittings presented here gets better 
results for exponential functions than any other, and even better 
to projection-corrected height than uncorrected. Actually, an ex¬ 
ponential ascent may be a characteristic of different ideal insta¬ 
bilities, e.g. "torus instability" at the first stages (Schrijver et al. 
2008; Gosain et al. 2012). Observing the filament expansion and 
its shape at several solar radii, it may reveal a kink instability, as 
in Torok & Kliem (2005). 

We checked the decay index n of the area from Sept 25 to 29, 
similar to the study by Zuccarello et al. (2014). The parameter 
Hcrit is always between 1.3 and 1.5 on the filament and emer¬ 
gence volume, which may indicate a torus instability. However, 
the lower part of the volume exhibits an important change of n in 
the area corresponding to the flux emergence, proving useful to 
identify ongoing instabilities. Importantly, Ucrit ranging the men¬ 
tioned values may also imply an exponential-to-linear expansion 
(Kliem & Torok 2006). 

After exploring the main models of eruption in the above 
paragraphs, ideal instabilities (torus instability and also kink) 
can be the most plausible mechanism. We may establish a possi¬ 
ble scenario as follows: the filament is destabilized, and a large 
magnetic flux emergence below the filament spine can help the 
destabilization, elevating against gravity and magnetic tension, 
increasing the magnetic pressure. 


The filament may be located close by a neutral point, which 
might be metastable (the paler shade of blue in Fig. 6, which is 
the minimum corresponding n in the Z direction). This point may 
allow the topology of the magnetic held to be maintained with¬ 
out 3-D reconnection with the lower part of the filament, and 
may help the torus instability to develop due to the flux emer¬ 
gence beneath this point. These kind of eruptions may be slow, 
leading to CMEs, and with one/several flux emergence patches 
underneath. 

4. Summary and conclusions 

This filament, located in the boundaries of the active region 
AR11850, evolves during almost five days of observations and 
eventually erupts. On Sept 29, 00:00 UT, a supergranular flux 
area started emerging and about 20 h later, the filament started 
rising. The most conspicuous rising phase lasted about 38 min 
and the Cl.2 flare onset started at 21:43 UT, while the two- 
ribbon arcade started developing around 15 min later. 

The height-time profile corresponds to an exponential pro¬ 
file, and the theoretical model that fits this smooth ascent most 
is the torus instability, plausibly related to a large magnetic flux 
emergence, and supported by the height-time profile and the suc¬ 
cession of events (first the fllament eruption and the flare after¬ 
wards). 
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